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Abstract 

Numerical studies of lattice-gas models are well suited to describe 
multi-adsorbate systems. One example is the underpotential deposi- 
tion of Cu on Au(lll) in the presence of sulfuric acid. Preliminary 
results from dynamic Monte Carlo simulations of the evolution of the 
adsorbed layer during potential-step experiments across phase transi- 
tions are presented for this particular system. The simulated current 
profiles reproduce a strong asymmetry seen in recent experiments. Ex- 
amination of the microscopic structures that occur during the simu- 
lated evolution processes raises questions that need to be investigated 
by further experimental and theoretical study. 

Keywords: Underpotential deposition. Multi-adsorbate systems. Lattice- 
gas model. Butler- Volmer approximation. Dynamic Monte Carlo sim- 
ulation. 



1 Introduction 



Underpotential deposition (UPD) provides a window of electrode potentials 
in which a monolayer or less of one metal can be deposited onto the surface 
of another metal. Recently, experiments that probe the dynamical prop- 
erties of the adsorption process have also added to our knowledge of 
UPD monolayers. Here we present preliminary results for dynamical simu- 
lations of a microscopic model which reproduce the results of potential-step 
experiments and make microscopic predictions that may be testable by mod- 
ern in-situ experimental techniques, such as scanning-probe microscopy and 
synchrotron X-ray scattering. 

The UPD of copper onto gold single-crystal electrodes in the presence of 
sulfate is a particularly well studied model system. See, for example, the re- 
view of the experimental literature in Ref. Q . For this system cyclic voltam- 
mograms, like the one shown in Fig. 1, show two current peaks corresponding 
to two different transitions in the surface structure |^]. Scanning-tunneling 
microscopy (STM) [Q, ^ and low-energy electron diffraction (LEED) [Q] 
studies indicate that these transitions occur between a full monolayer (ML) 
of Cu at more negative potentials, an ordered (\/3 x \/3) mixed copper 
and sulfate phase at intermediate potentials, and a disordered low-coverage 
phase at more positive potentials. In situ X-ray scattering confirms that 
the ordered phase consists of a 2/3 ML of copper in a honeycomb pattern, 
with a 1/3 ML of sulfate occupying the remaining adsorption sites [P, as 
originally proposed by Huckaby and Blum [^. 

An effective theoretical approach for understanding a multi-species layer 
chemisorbed on a surface is to study a statistical-mechanical lattice-gas 
model numerically. Such models describe the adsorbate layer by a Hamil- 
tonian which gives the energies of different adsorbate configurations. Early 
studies [0, P together with later applications to problems such as the 
adsorption of urea on Pt(lOO) [10, 11 1 and the UPD of Cu with sulfate on 
Au(lll) @, 0], in or near equilibrium, are discussed in recent review articles 
[ll2| , [l^ , |l4| . In particular, that work provided estimates for effective lateral 
interaction energies, based on comparison of the theoretical predictions with 
experimentally observed equilibrium adsorption isotherms and voltammetric 
currents for small potential sweep rates. 

Dynamical properties of the copper UPD system have recently been in- 
vestigated by Holzle, Retter, and Kolb (HRK) They measured current 
transients in potential-step experiments performed at both transitions and 
with both positive-going and negative-going steps. In this article we re- 
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port current transients and time-dependent microscopic adlayer structures 
observed in dynamic Monte Carlo simulations of a model UPD system un- 
der conditions intended to reproduce those experiments. In a previous 
study Ijl^, our model agreed qualitatively with the experiments only for 
steps across one transition, the mixed-disordered phase transition marked A 
and A' in Fig. 1. Reproducing the behavior seen for potential steps across 
both transitions represents a significant refinement of the lattice gas model. 
Our microscopic predictions indicate the desirability of time-resolved, struc- 
turally sensitive, in situ experimental studies of this system. 



2 Lattice-gas Model 

Our theoretical lattice-gas model treats copper and sulfate as interacting 
particles that compete for the same adsorption sites. For the Au(lll) surface 
this is a triangular lattice of threefold hollow sites Q . The model is defined 
by the effective lattice-gas Hamiltonian, 
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Here equals unity when lattice site i is occupied by species X, with 
X = C for copper and X = S for sulfate; otherwise equals zero. This 
third state represents sites solvated by water. The effective interaction for 
n-th neighbor adsorbates X and Y is ^I'xY' ^^'^ ^fe^) ^^P^^sents the sum 
over all n-th order neighbor pairs. We also include a three-body interaction, 
$gg, for equilateral triangles of sulfate adsorbed on second-neighbor sites. 
The sum over all possible triangles is denoted by J2a- Finally, fix is the 
electrochemical potential for X, and J2i is the sum over all lattice sites. 
Attractive interactions are represented by positive quantities, and positive 
electrochemical potentials indicate a preference for adsorption in the absence 
of lateral interactions. 

(n) 

The quantities $XY lattice-gas model are effective parameters rep- 

resenting the joint effect of competing, unspecified interaction potentials. 
Ground-state calculations can be used to determine multi-dimensional re- 
gions of parameter space that give phase diagrams that agree with the phases 
observed experimentally. Fitting the model to a particular system involves 
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an iterative process of selecting experiments and adjusting the parameters 
to achieve improved agreement. The expectation is that the model will more 
accurately represent the real chemical system as it correctly reproduces and 
predicts an ever wider range of experimental observations. Detailed discus- 
sion of the methods used to estimate the effective interaction constants are 



given in Refs. |2||ll|, |12|. 



The effective interactions used in the present study are illustrated in 
Fig. 2. There, and elsewhere in this paper, unoccupied adsorption sites are 
represented by an open circle (o), sites occupied by copper are represented 
by a filled circle (•), and those occupied by sulfate are represented by a 
triangle (A). These interactions have been adapted from those of Ref. [^. 
Changes in the sulfate-sulfate interactions removed from the ground-state di- 
agram three experimentally unobserved low-temperature phases which were 
overlooked in Ref. [||. Strengthening the attractive next-nearest neighbor 
copper-copper interaction, ^qq, was required to make the phase transition 
between the copper monolayer and the mixed phase (peaks B and B' in 
Fig. 1) first order and produce current transients for this transition in qual- 
itative agreement with experiments (described in Sec. 4). 

The ground-state diagram for the interactions of Fig. 2 is presented in 
Fig. 3, and the geometries of the main observed phases are illustrated in 
Fig. 4. Here the phases are labeled by the length of their basis vectors and 
subscripted (superscripted) by the fraction of sites filled by copper (sulfate) . 

The electrochemical potential of species X is controlled through its bulk 
solution activity and the electrode potential, E. In the weak-solution limit, 
the activity equals the solute molar concentration [X]. Then the electro- 
chemical potential, fix, is 

[X] 



fix = fix + RTln^ - zxFE, 



where zx is the effective electrovalence of X, F is Faraday's constant, R is 
the molar gas constant, T is temperature, and the reference state is /i^) [X]*^, 
and = 0. For fixed [S] and [C], Eq. (|2|) for fis and fie form the parametric 
representation of a line in the ground-state diagram. The electrode potential 
E is the parameter that determines the location along this isotherm, whose 
slope is zs/zc- We have chosen an isotherm such that the potential difference 
between the two transitions is the same as observed by HRK. Its location is 
given by the dotted line in Fig. 3 for < E < 350mV vs Ag/AgCl. 
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3 Dynamic Monte Carlo Method 



The equilibrium properties of a lattice-gas Hamiltonian can often only be 
practically evaluated using Monte Carlo methods, by which thermodynamic 
average values are calculated from a representative sampling of the equilib- 
rium states. This sampling is generated using the principle of detailed bal- 
ance: in equilibrium the rate of going from configuration A to configuration 
B is equal to the rate of going from B to A. Many Monte Carlo techniques 
generate new configurations through a series of Monte Carlo steps, each step 
consisting of a proposed small change in the configuration which is accepted 
or rejected according to an acceptance probability that must satisfy detailed 
balance. 

A lattice-gas Hamiltonian merely provides the free energies correspond- 
ing to different geometric configurations of the adsorbate layer. In contrast 
to a quantum-mechanical Hamiltonian, it does not possess an intrinsic dy- 
namic. However, when Monte Carlo steps are associated with a discrete unit 
of time, the simulation can be viewed as a "movie" of a dynamic process. 

Many different Monte Carlo algorithms can be designed to describe the 
equilibrium properties of a particular system. In general such an algorithm 
does not realistically describe the dynamics of a particular chemical sys- 
tem. But if care is taken to include only dynamically realistic moves in the 
algorithm, it is possible to construct a Monte Carlo dynamic which simu- 
lates that of the real system under consideration. This should be the case 
when the real dynamics is well approximated by a stochastic process on the 
timescales of interest. The example relevant to our present discussion is 
the thermally activated motion of particles adsorbing onto, desorbing from, 
and diffusing between adsorption sites on the electrode surface. This view 
of Monte Carlo simulation can be combined with the lattice-gas Hamilto- 
nian to construct a stochastic dynamic describing the individual motions 
of specific adsorbate particles. An introduction to Monte Carlo methods 
for equilibrium and nonequilibrium problems in surface electrochemistry is 
given in Ref. []T6|| . 

Here we construct a microscopic dynamic which includes adsorption, 
desorption, and one-step lateral diffusion of both adsorbate species. The 
rate of each process is determined using a simple Arrhenius picture, with 
all processes having the same pre-exponential factor vq. The free energy 
associated with an adsorbed particle of species X occupying site i for a 
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specific configuration of neighbors is given by 



W,,(X;r) = -^ 



(n) (n) 



'^x,s4s E c'c' - /^x , (3) 

A(i) 



where Sj"^) represents the sum over all adsorption sites j that are n-th 
neighbors of site i, and J2A{i) is the sum over all second-neighbor equilateral 
triangles involving site i. The factor 5x,s is unity when X = S and zero 
otherwise. The index T runs over all possible arrangements of neighboring 
adsorbate particles within the maximum interaction range from site i. When 
the site is empty and the ion is in solution, the free energy is defined to be 
zero. Using X = to denote an unoccupied site, TCi{0;T) = 0, regardless of 
the arrangement of the neighbors. 

The energy barriers are assumed to vary according to the Butler- Volmer 



approximation |17]. Specifically, consider the adsorption-desorption pro- 
cesses for a particular site. A schematic free-energy surface is shown in 
Fig. 5(a), assuming three different arrangements of neighbors, F. If 'Hj(X; F) = 
0, there is no preference toward adsorption, but there is still a transition- 
state free energy, Aq(X). The transition-state free energy for this case is 
used as the reference value for the adsorption transition energy for all con- 
figurations of the neighbors. If a site has no neighbors within the maximum 
interaction range, the free energy of the adsorbed particle is controlled by the 
electrochemical potential and equals —fix- In the Butler- Volmer approxi- 
mation, the transition energy for adsorption to that site is then lowered 
from Aq(X) by an amount a fix- Here a is the symmetry factor, which is 
often close to 1/2. For the present model, lateral interactions with neigh- 
boring ions are also important. If the net lateral interactions are attractive, 
the free energy corresponding to the adsorbed state lies below —p-x, as il- 
lustrated by the lowest adsorbed free-energy level in Fig. 5(a). If the net 
lateral interactions are repulsive, the free energy is raised above —fix- Our 
generalization of the Butler- Volmer approximation uses the total free energy 
gained to calculate the free-energy barrier for adsorption, specifically 

A*(X;F) = AS(X)+aW,(X;F). (4) 

Using the same transition-state energy A*(X; F) for the adsorption and des- 
orption process is essential for Monte Carlo simulations in order to satisfy 
detailed balance and thus ensure eventual approach to thermodynamic equi- 
librium. In the model developed here, an ion approaching the surface from 
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the solution sees different free-energy barriers for sites with different arrange- 
ments of neighbors. Adsorption on a site is faster when its neighbors help 
lower the free energy of the adsorbed state; conversely, desorption from such 
a site is relatively slow. For multi-component adsorption, the adsorption of 
one species may be much faster than that of the other at a particular site. 

For lateral diffusion from one adsorption site to another, a slightly mod- 
ified approach is needed. The situation for two neighboring sites, both 
with favorable adsorption free energies, is shown in Fig. 5(b). The diffi- 
cultly is that detailed balance requires a single transition-state free energy, 
A*j(X;r), between the two sites. An obvious solution would be to adjust 
the energy barrier by the average of the free energies for the two sites, 
{ni{X;T) + nj{X;T))/2 |§. Unfortunately this could lead to Wi(X;r) < 
A*j(X;T) < T-Cj(X;T), which violates the usual definition of a transition- 
state free energy. To avoid this, we use the least favorable adsorption free 
energy for calculating the transition-state free energy for diffusion between 
nearest-neighbor sites i and j: 

A*^.(X;r) = AS(X) + iMax(W,(X;r),W,(X;r)) . (5) 

Violations of the transition-state picture are minimized because they in- 
volve very unfavorable sites, which are only rarely occupied. The barrier 
for diffusion to each neighbor is calculated individually; diffusion into a less 
favorable site will be relatively slow. Diffusion to an occupied site is for- 
bidden and has a rate identically zero. For our simulations, we choose the 
transition-state free energies as Aq(C) = 50kJ/mol for copper adsorption, 
A^(S) = 55kJ/mol for sulfate adsorption, and A5(C) = A^(S) = 35kJ/mol 
for diffusion of both species. These choices make adsorption faster for copper 
than sulfate and allow relatively long diffusion paths across the surface since 
diffusion is faster than adsorption. We take the symmetry factor a = 1/2. 
The rates described above can be used to implement a rejection- free 



Monte Carlo algorithm |16, A list of all the single-ion processes is 

maintained. Each process is weighted by its rate, so that making a random 
choice from the list is biased towards the faster processes. Each simulation 
step consists of choosing a process from the weighted list, making the ap- 
propriate change to the system, and recalculating the weights for the list. 
At each step the simulation clock is updated by a random number that rep- 
resents the elapsed time since the last move. For a given configuration the 
overall rate for the system to change by some microscopic process is the sum 
of all the individual microscopic process rates. Since the decay of a particu- 
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lar state has a constant rate, the elapsed time for the actual change has an 
exponential distribution. Results for a simplified model were presented in 
Ref. 

Our simulations were conducted on 30 x 30 parallelograms of a triangular 
lattice with periodic boundary conditions. A particular potential-step simu- 
lation involves equilibration for a fixed time at one potential. Then, at time 
t = 0, the applied potential is changed to its final value, and measurements 
of the partial coverage due to each species are made at fixed simulation times 
after the quench. Results reported here are averaged over at least 1000 in- 
dividual simulations. Currents are calculated by numerical differentiation 
of the coverage with respect to time, with electrovalences of zs = —0.9 and 
zc = +1.54 Q. The natural simulation time unit is Monte Carlo steps per 
site (MCSS), but choosing uq = 2 x 10^ MCSSs""^ makes the scale of simu- 
lated currents agree roughly with experiments. The rejection- free algorithm 
is essential, since achieving a time scale of 10^ MCSS would be prohibitive 
for a traditional Monte Carlo algorithm. Even with the rejection-free algo- 
rithm, some of the individual experiments took 15 minutes to simulate on a 
DEC- Alpha workstation. 

4 Results 

The agreement between our simulations and the experiments of HRK |||] is 
encouraging (see Figs. 6 and 7). The current profiles, consisting of a rapid 
initial decay followed by a broad maximum at a later time, are character- 
istic of the decay of an adsorbed layer that is metastable after a potential 
step across a first-order phase transition. The initial current transient cor- 
responds to the rapid relaxation to a long-lived metastable state which has 
only a small concentration difference from the adsorbate layer as it existed 
before the potential was changed. The second, broader current maximum 
occurs as domains of the equilibrium phase nucleate and grow. This behavior 
is seen for three of the four potential steps. 

The simulated current transients for potential steps around the transition 
between the (1 x 1)5 and (\/3 x \/3)2/3 phases, labelled B and B' in Figs. 1 
and 3, are shown in Fig. 6. Both reflect one-step nucleation and growth 
mechanisms. Following the trend seen in the experimental data, the broad 
maximum becomes stronger and the maximum shifts to earlier times as the 
size of the step is increased. However, the time corresponding to the second 
maximum appears to depend more strongly on the step size than observed 
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in the experiments. 

The current transients observed for the transition between the (\/3 x 
\/3)2/3 and the disordered fow-coverage phases, labelled A and A' in Figs. 1 
and 3, are very different for positive and negative potential steps, as seen 
in Fig. 7. This was also observed by HRK. When the potential is increased 
from the ordered-phase region into the disordered-phase region, a step corre- 
sponding to A', the current profile has the broad maximum, similar to those 
shown in Fig. 6. Again, the qualitative trend seen for the strength and loca- 
tion of the second current maximum as the size of the step is changed agrees 
with the experiments. In contrast to our previous simulations, which were 
performed for a system with different interactions [15|, this current is caused 
by the decay of the metastable (\/3x \/3)2/3 phase to the (\/3x -v/3)j^^2 phase, 
which is also metastable. The decay of this second metastable phase to the 
final low-coverage phase takes place at much later times, but causes a cur- 
rent transient with a much smaller amplitude because the copper and sulfate 
currents have opposite sign. Both metastable phases decay by a nucleation 
and growth mechanism. 

The current profiles observed for potential steps in the negative direction, 
corresponding to peak A, are qualitatively different than the other profiles. 
They are monotonically decreasing and do not correspond to a simple func- 
tional form. A sequence of configurations after a potential step to 20 mV 
below the transition is shown in Fig. 8. Because the applied potentials of 
HRK have been used, the equilibration occurs in a disordered phase with 
a relatively high sulfate coverage of approximately 0.25 ML. After the po- 
tential step, copper adsorbs quickly onto the surface. The sulfate already 
present collapses with this copper to form a domain of the (\/3 x \/3)J/3 

phase. Later this fills in with copper, forming the equilibrium (\/3 x \/3)2/3 
phase. The adsorption of sulfate is an even slower process, which controls the 
shrinking of the domains of unoccupied sites. This process is different from 
our previous simulations , where no domain formation was observed. In- 
stead, copper adsorbed into a loose network with open sites favorable for 
sulfate adsorption. 

The fact that different dynamical paths can give current profiles quali- 
tatively similar to those observed experimentally indicates that more study 
of the copper UPD system is desirable. Particularly useful would be time- 
resolved experiments that can probe the microscopic structure of the ad- 
sorbed layer. For the simulations, the factors controlling the different dy- 
namical behaviors are not clearly understood and a more thorough investi- 
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gation is needed. 



5 Conclusions 

Preliminary dynamic Monte Carlo simulation results are reported for po- 
tential step experiments for UPD of copper with sulfate on Au(lll). The 
simulated current transients display the same general characteristics as those 
observed experimentally by Holzle, Retter, and Kolb Q. The current profiles 
for potential steps across the transition separating the copper monolayer and 
the ordered mixed phase (the transition marked B and B' in Fig. 3) show 
a sharp initial transient followed by a broad second maximum that shifts 
to earlier times and larger amplitudes as the size of the step is increased. 
This contrasts with the asymmetry of the current profiles for potential steps 
across the transition separating the ordered mixed phase from the disor- 
dered low-coverage phase (the transition marked A and A' in Fig. 3), where 
the second maximum is only seen for positive- going steps. Negative poten- 
tial steps produce a monotonically decreasing current. We find that repro- 
ducing the profiles of HRK is not sufficient for determining the underlying 
dynamical process. For two slightly different models we have observed dif- 
ferent processes occurring after the step. For both the negative-going and 
positive-going steps at this transition (A and A'), an intermediate ordered 
phase with only 1/3 monolayer of copper is observed in the model presented 
here that was not observed for an earlier model. In addition, domain for- 
mation is seen for negative steps. This was not observed in our simulations 
of the earlier model. Experiments which can dynamically monitor the mi- 
croscopic structure of the adsorbate layer during the relaxation after the 
potential step would be useful in guiding further model development. 
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Electrode Potential (mV vs. Ag/AgCl) 



Figure 1: Cyclic voltammogram for Cu UPD onto Au(lll) surface, after 
Ref. Q. The two peaks correspond to two first-order phase transitions in 
the adsorbed layer. The arrows indicate the direction of the potential scan. 
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Figure 2: Diagrams illustrating the relative adsorbate positions correspond- 

(n) 

ing to the n-th neighbor interactions, <&XY' ™ lattice-gas model. Here Cu 
is represented by a filled circle (•), sulfate is represented by a triangle (A), 
and unoccupied sites are open circles (o). The interaction energy in kJ/mol 
is given below each diagram. Positive values indicate attractive interactions. 
First-neighbor sulfates are subject to an infinitely strong repulsion, i.e. con- 
figurations with first-neighbor sulfates are forbidden. The lower right-most 
diagram illustrates the three-body sulfate interaction. Note that some of 
the interactions are different from those used in Ref. 
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Figure 3: The ground-state phase diagram (sohd lines) for the lattice-gas 
model with interactions as given in Fig. 2. The axes represent the elec- 
trochemical potentials of the sulfate and copper, fls and /Uc, respectively. 
The isotherm (dotted line) used for the simulated potential-step experiments 
was chosen to match the transition potentials observed in the experiments 
of HRK jl| . The geometries of the different phases are illustrated in Fig. 4. 
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Figure 4: Important phases observed for the present lattice-gas model of 
copper and sulfate adsorption on Au(lll). Cu is represented as a filled circle 
(•), sulfate as a triangle (A), and empty sites as an unfilled circle (o). The 
label indicates the length of the unit vectors. The subscript and superscript 
represent the fraction of sites occupied by copper and sulfate, respectively. 
Two sulfate-free copper phases, (-\/3 x -\/3)5/3 and (\/3 x \/3)2/3) occur for 
this model at large negative /y,s, but arc not shown. Their geometries are 
easily deduced from the phases that are shown. 
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Figure 5: Butler- Volmer type free-energy barrier scheme. See text for de- 
tails. 




Figure 6: Simulated current transients after potential steps across the tran- 
sition between the (1 x 1)5 and (\/3 x \/3)2/3 phases. Positive (top) and neg- 
ative (bottom) potential steps, B' and B, respectively, of HRK Q are shown. 
The shape of the current profile corresponds to the decay of a metastable 
adsorbed layer. The trend for the broad maxima to shift to earlier times and 
larger amplitudes as the potential step is increased agrees with the experi- 
ments. The legend indicates the size of the potential step past the transition, 
l-Efinai — -^'transition I • This figure should be compared with Fig. 3 of Ref. |||]. 




Figure 7: Simulated current transients after potential steps across the 
disorder-mixed transition for positive (top) and negative (bottom) potential 
steps, A' and A, respectively, of HRK The strong asymmetry in cur- 
rent profiles for the two step directions is seen experimentally. The broad 
maximum is characteristic of the decay of a metastable adsorbed layer, the 
monotonically decreasing initial current reflects the unstable nature of the 
low-coverage surface after a negative potential step. The legend indicates 
the size of the potential step past the transition, l-Egnai — -^'transition I • This 
figure should be compared with Fig. 3 of Ref. 






Figure 8: A series of snapshots after a negative-going potential step to 20mV 
below the transition between the low coverage and mixed layers. After the 
step, a fraction of the sulfate desorbs, but the remaining sulfate combines 
with the newly adsorbed copper to form a loose domain. With time this 
domain fills in and grows. The current transient is monotonically decreasing, 
as shown in Fig. 7, A. Movies of the dynamics after the quench are available 
in MPEG format at http : //www, scri . f su. edu/~rikvold| 



